Diffusion-reaction in thermal growth of silicon oxide films on Si 
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The thermal growth of silicon oxide films on Si in dry O2 is modelled as a dynamical system, 
assuming that it is basically a diffusion-reaction phenomenon. Relevant findings of the last decade 
are incorporated, as structure and composition of the oxide/Si interface and O2 transport and 
reaction at initial stages of growth. The present model departs from the well established Deal and 
Grove framework (Deal, B.E. and Grove, A. S. General Relationship for the Thermal Oxidation 
of Silicon J. Appl. Phys. 36, 3770-3778 (1965)) indicating that its basic assumptions, steady-state 
regime and reaction between O2 and Si at a sharp oxide/Si interface are only attained asymptotically. 
Experimental growth kinetics by various authors, obtained for a wide range of growth parameters 
are shown to collapse into one single curve when the scaling properties of this model equations are 
explored. 



Silicon oxide films thermally grown on single- 
crystalline Si (c-Si) wafers in dry O2 are still the most 
common materials used as gate dielectrics in Si-based 
metal-oxide-semiconductor (MOS) structures. As device 
dimensions shrink below 0.25 /im, the oxide films are 
forced to scale down to thicknesses of 5 nm and less. Re- 
liability of highly integrated Si devices is critically depen- 
dent on the characteristics of the vitreous oxide film, like 
thickness uniformity, defect density, dielectric strength, 
and others, as well as on those of the oxide/Si interface, 
like roughness, electronic states density and others 

The dielectric performance of ultrathin films of silicon 
oxide thermally grown on Si and the structural and elec- 
tronic properties of the oxide/Si interface have been in- 
tensively studied experimentally, theoretically, and com- 
putationally in synergism with the development of 
the semiconductor industry. The understanding of film 
growth kinetics, however, did not make any significant 
progress since the early model of Deal and Grove ||, 
which subsisted as the only fundamental approach. The 
Deal and Grove model, and more specifically the result- 
ing linear-parabolic law, agrees with the observed growth 
kinetics in dry O2 only above a certain thickness HI). 
This is an expected fact from the model itself, since its 
basic assumption is that growth is promoted by inter- 
stitial diffusion of the oxidant species, the O2 molecule, 
through a previously grown oxide layer, thick enough to 
guarantee steady state regime. Besides, reaction between 
O2 and Si is assumed to take place solely at an abrupt 
oxide/Si interface, producing stoichiometric Si02- This 
initial oxide thickness was estimated from experi- 
mental data as being between 20 and 30 nm, well above 
the thickness range of interest for present and future use 
in microelectronics Si-devices. Description of the growth 
kinetics in the lower thickness range was addressed by 
many authors, within the Deal and Grove model frame- 
work, adding new terms to the linear-parabolic expres- 
sion in order to fit experimental data Although 
fitting the experimental curves, and so providing useful 



analytical expressions capable of reproducing the whole 
thickness interval, the extra terms added to the Deal and 
Grove expression did not have a well-defined physical 
meaning, even though their dependence on some process- 
ing parameters, like temperature for instance, could be 
explored || . None of these models has been shown to be 
clearly correct and none of them has gained widespread 
acceptance. One is left then with several empirical ex- 
pressions which can be used to model the growth kinetics 
in dry O2 in the thin and ultrathin film regimes below 
20 nm, as well as in further stages of growth, but with 
no satisfactory physical explanation H. According to 
B. Deal |l(| "such an expression may well fit the data, 
but provides little guidance or insight into mechanisms 
involved or effects of process variations" . 

Recent investigations showed that as film thickness de- 
creases below 20 nm, the contribution to film growth 
due to reaction away from the interface region becomes 
increasingly significant |ll|-|l6[]. Furthermore, there are 
strong theoretical and experimental evidences of the exis- 
tence of a reactive layer formed by sub-oxides (also called 
Si-excess) near the oxide/Si interface [T^-pOl]. 

In the present work the thermal growth of silicon oxide 
films on Si is modelled as a dynamical system, assum- 
ing that it is basically a diffusion-reaction phenomenon. 
Since steady-state regime is not imposed, an initial ox- 
ide thickness is not required. Therefore, the model is 
expected to describe the whole oxide thickness interval. 
Diffusion-reaction equations have been used to describe 
different systems [^l]-^3| but, to our knowledge, this is 
the first time they are applied to silicon oxide growth 
kinetics. The diffusing species is taken to be O2, to 
model what has been largely demonstrated by isotopic 
substitution experiments |ll|,|l2|,^4|-|2^] . Growth is pro- 
moted by reaction of O2 with Si, not necessarily at the 
oxide/Si interface, but wherever the two species meet. 
The growth kinetics can then be obtained at any tem- 
perature by specifying the diffusivity of the O2 molecule 
in the silica network, D, the reaction rate between O2 and 
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Si, k, and the O2 pressure in the gas phase. Growth in 
one dimension is considered, which is the most common 
experimental situation, meaning that as a face of c-Si 
is exposed to O2, the silica film grows in the direction 
perpendicular to this face. The proposed description for 
the growth kinetics is contained in the following coupled 
partial differential equations: 
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where pi — c^/cSf* is the concentration of the i-species 
(i = Si, O2) in the solid phase, in units of number of 
atoms per unit volume, c$, normalised by the concentra- 
tion of Si in c-Si, c b gf k . Thus, the state of any species 
at any time is characterised by a density function p(x, t), 
where x is the coordinate in the direction of growth, x = 
indicating the surface of the sample, and t is the elapsed 
growth time. The first equation describes the O2 concen- 
tration rate and it is essentially a balance equation. At a 
given position the concentration of oxygen may vary due 
to the net flux of O2 through fully (SiC>2 in oxide bulk) or 
partially (SiO^, x < 2, near the interface fl jl^ , pj| , p0| ) ox- 
idised Si, modelled by the diffusion term with a diffusion 
constant D, or due to reaction to produce silicon oxide, 
characterised by a reaction rate constant k. The sec- 
ond equation describes the Si concentration rate, which 
is only due to reaction with O2, implying that Si is an 
immobile species p4| . The silicon oxide concentration, 
Poxide, does not explicitly appear in these equations, but 
it is in direct connection with psi through the expression 



Poxide(x, t) = 1 - pSi(x, t) 



(2) 



reflecting conservation of Si species in the Si + O2 — > 
SiC>2 reaction. It is assumed that the O2 diffusivity at 
a certain temperature does not change during oxidation, 
and therefore D is constant throughout the whole growth 
process. Furthermore, based on the fact that oxygen dif- 
fuses mainly through oxide, one neglects volume changes 
due to the above chemical reaction. In order to complete 
the description in terms of equations ([!]), it is necessary 
to provide the initial and boundary conditions. They are 



p Si (x,0) = l Vx>0 



(3) 



Po 2 (0,t) = 



Cgasfv 
„bulk 



Vi > 



where c gas is the O2 concentration in the gas phase and 
f v is the ratio between the accessible free volume for O2 
in the silica network and the unit volume of the solid || . 
Initially there is a pure, c-Si substrate from x = to 00, 
whose surface is exposed to a gaseous medium, O2, with 
a relative concentration at this surface represented in the 



second condition by 



Reaction of these two species 



produces silicon oxide, not necessarily stoichiometric, as 



the solid-phase reaction Si + O2 — > SiC>2 may occur in 
several steps. As time proceeds, O2 diffuses through the 
oxide network, reaches Si, reacts producing more oxide, 
and consequently pushes the oxide / Si interface further in 
the positive x region. Different values of D, k and c gas 
imply different kinetics. However, the values of D and 
k may be used to define "natural" units for the system, 
resulting in a set of adimensional equations with only one 
parameter, namely c gas p3[ . Under the following trans- 
formations: 
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<pi(u,r) = pi(x,t) for i = Si, O2, oxide 
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the resulting adimensional set of differential equations 
and initial and boundary conditions are 
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leaving 



■gas 



as the only parameter. This is a powerful 



statement. It predicts that in all experiments performed 
at the same gas pressure, oxide growth kinetics and O2 
and Si concentrations versus depth curves (profiles) at 
different temperatures collapse to the same curves, pro- 
vided that appropriate units are used: 1/k for time and 
y/D/k for thickness. 

In order to obtain profiles and growth kinetics, eqs.((^) 
must be solved. Those are partial, nonlinear, coupled 
differential equations and cannot be analytically solved 
at all times. On the other hand, they can be easily solved 
by finite differences methods, iterating for as long as it 
takes to reproduce the growth time of interest. The cal- 
culation gives then the profiles <j>i(u, t) of all species and 
the kinetics. Typical profiles at different growth times 
for a given value of c gas are shown in Figure 1. The O2 
profiles are not straight lines as they are assumed to be 
in the steady-state regime || . Strictly speaking, steady- 
state is never attained as long as the oxide/Si interface 
moves deeper into the c-Si substrate. Furthermore, this 
interface is not sharp, in accordance with a wealth of 
experimental evidences [ p0|p5|j2^ |. Still in strict accor- 
dance with theoretical predictions ]4| JT7|Jl^| and experi- 
mental observations ]l9| , p0| , f25| -27|, a finite width region 
formed by Si in various lower oxidation states (+1, +2, 
+3, sub-oxides or Si-excess) early develops between the 
stoichiometric SiC>2 film and the c-Si substrate. In this 
way, 02-Si reaction can take place in the whole region 
where non-fully oxidised Si is available. 
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(i.e. different D and k) when natural units of the sys- 
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FIG. 1. Calculated Si and O2 profiles in the solid phase 
at various temperatures and one pressure, for different oxida- 
tion times. The solid, dashed, and dot-dashed lines represent, 
respectively, the profiles at increasing oxidation times. 

Fig. 1 shows that the two hypotheses of the Deal 
and Grove model are not valid: i) Steady-state regime 
is never attained; ii) reaction between O2 and Si leading 
to growth does not take place at a sharp oxide/Si inter- 
face. For high enough temperatures and pressures, and 
long enough times, the O2 profile approaches a straight 
line, and the thickness of the near-interface, sub-oxide 
region is much smaller than the oxide thickness. This 
is the observed situation at oxide thicknesses well above 
20 nm. In this limit, steady state and abrupt interface 
are plausible approximations. In other words, asymptot- 
ically the Deal and Grove assumptions are valid, whereas 
at initial stages they are not. The kinetics in the initial 
growth regime, which was described as "anomalous" in 
the Deal and Grove framework, naturally emerges from 
the present model as follows. Let us consider the quan- 
tity that is mainly focused in practical situations: the 
oxide thickness x- We calculate it through the relation: 



X(r) = 



4>oxide{u,T)du= I (l—(f>Si(u,T))du (6) 
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The oxide growth kinetics x( T ) is completely determined 
by <t>Si{ u , T ), whose evolution is ruled by eqs.(||) and 
therefore is determined by a unique parameter c gas , the 
concentration of O2 in the gas phase. However, due to 
the particular form of the differential equations, together 
with boundary and initial conditions, the system presents 
a further scaling property. Suppose two different growth 
kinetics, with O2 concentrations in the gas phase given 

by c^as and cfas 
apply 



then the following relations 
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tern are used, kinetic curves for different 



collapse 



to one single curve due to the symmetry present in the 
form of the diffusion-reaction model equations together 
with initial and boundary conditions. In summary, the 
theoretical prediction is that all kinetic curves reduce to 
a single one. 

These predictions were tested with kinetic curves taken 
from references f|p7|-p9f , corresponding to oxidations in 
Joule-effect heated and rapid thermal processing furnaces 
in a variety of growth parameters (temperature, O2 pres- 
sure, and time) leading to oxide thicknesses ranging from 
1 to 100 nm. In Fig. 2 we show the collapse of the exper- 
imental kinetics (symbols, colours) into the theoretical 
curve (line, black). For the kinetics curves from reference 
p8fl , with five points in each curve, we first used the scal- 
ing predicted by eqs.(Q), obtaining four different curves, 
one for each temperature. We then scaled these curves 
using natural units for x and t according to eqs. (@) until 
they overlapped the theoretical curve. The curves from 
reference || contained much more experimental points, 
extending from 0.1 to 1000 min, and from 1 to 80 nm. 
For these curves the pressure scaling predicted by eqs.(0) 
is not perfect at long growth times, probably due to an 
increase in diffusivity. They were then collapsed one by 
one into the theoretical curve using eqs.(Q) first, and then 
eqs.(0). The same procedure was used to collapse the ki- 
netics from references |27|,^1 into the theoretical curve. 
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FIG. 2. Calculated growth kinetics (solid line, black) and 
experimental growth kinetics from the indicated references 
(symbols, colours) collapsed into the calculated curve. 

The asymptotic depart from the theoretical curve is 
an expected result, since the present model assumes a 
constant diffusivity throughout the whole growth pro- 
cess, despite many evidences that this is not the case: i) 
transport of the diffusing species through the O-excess 
region near the surface [l2jJl~^ , ^G| . |30] | , thr ough the Si - 



excess region near the oxide/ Si interface 15 18j,|2 



Therefore, besides collapsing for different temperatures and through the stoichiometric Si02 in the bulk of the 
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growing oxide film have different diffusivities; ii) the 
growth of the defective, near-surface and near-interface 
regions saturate within a few nm |l^,|f^,^,|^,|o) , while 
the bulk, stoichiometric oxide grows continuously. So, 
interstitial diffusion of O2 through the bulk of the grow- 
ing oxide gives a relative contribution that increases as 
the width of this region becomes dominant; iii) oxida- 
tions performed in dry O2 flow during long time inter- 
vals (hours) may suffer from water vapor contamination, 
which accelerates the oxide growth [§[l3| . 

Taking the dimensional oxide thickness X(i) = 

X( r )\/f> and dennin g Po = Po a {0,t) = the 
asymptotic solution to eqs.(ra) predicts the inverse growth 
rate dt/dX = jj^X, with D constant, which is very 
similar to the asymptotic approximation of the Deal and 
Grove model § dt/dX = ^X. In Fig. 3 the inverse 
growth rate is plotted as a function of oxide film thickness 
for different experimental kinetics ||. 



to a quantitative description of the complete growth ki- 
netics. 

The present approach opens the possibility of incor- 
porating new and determinant facts on O2 transport at 
initial stages of thermal oxide growth and on structure 
of the oxide/Si interface, that were discovered after the 
proposal of Deal and Grove || and not yet integrated 
into growth kinetics models. Some were discussed here, 
like the graded, Si-excess nature of the oxide/Si interface 
and the incorporation of freshly arriving O2 into different 
regions of the growing oxide, while several others may 
also be included in a diffusion-reaction approach, such 
as the role played by oxygen excess centers j30) in the 
near-surface region on oxygen transport which has been 
recently investigated [ ^6| . 




X(nm) 

FIG. 3. Experimental inverse growth rates [6] and theoret- 
ical inverse growth rate (inset). 

A deviation from a straight line in the thick film range 
is apparent, corresponding to a variable, increasing dif- 
fusivity, contrary to the theoretical prediction shown in 
the inset. 

Thus, the thermal growth of silicon oxide films on c- 
Si in dry O2 is promoted by diffusion of O2 through the 
growing oxide with variable diffusivity. Similar physi- 
cal arguments could support a variable diffusion rate. 
Nevertheless, even using constant, effective values for the 
diffusivity and reaction rate, modelling the growth as a 
diffusion-reaction process where conditions such as initial 
thickness, steady-state regime, and reaction solely at an 
abrupt interface are abandoned allows to reproduce the 
observed kinetics, especially in the lower thickness (ultra- 
thin film) regime which is the one of present and future 
practical interest. Realistic model evaluations of variable 
diffusivity and reaction rate stay beyond the scope of the 
present work, although this seems to be the only route 
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